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ABSTRACT 


A representative sample of lunar rilles observed on Lunar Orbiter IV and 
V photographs is presented. A method of classification of rilles is described 
which will provide a means of quantitatively grouping rilles of similar plane- 
metric shape. The method is based on a finite Fourier approximation of the 
edge of the rille. It provides an adequate approximation to all planemetric 
characteristics of the rille except rille recurvature. 


INTRODUCTION 


Lunar rilles have been subjectively classified as sinuous, normal, 
straight, arcuate, and branching. While there seems to be general agreement 
for the origin of many rilles, the origin of sinuous rilles has always been a 
subject of considerable controversy. The diversity of theories proposed for 
the production of sinuous rilles includes: surface and subsurface erosion by 

flowing water, erosion by ash flows, lava channel formation, lava tube col- 
lapse, magma intrusion followed by collapse, gas venting along subsurface 
fractures and intersection of fracture patterns . The great diversity in 
theories may indicate real differences in shape and origin of different sinu- 
ous rilles. If so, the subjective definition of sinuous rilles is much too 
general. Alternatively, all sinuous rilles may be similar in shape which 
might favor a single theory for all sinuous rilles. It is clear, therefore, 
that an objective classification of lunar rilles in general and sinuous rilles 
in particular is required if the origin of rilles is to be studied on a 
quantitative basis. 

The purpose of this paper is to present a catalog of a representative, 
group of lunar rilles and a brief description of a quantitative method of 
rille classification. The catalog will be used to obtain a subjective classi- 
fication of lunar rilles from workers actively engaged in lunar studies. The 
subjective classification can then be compared to the results of an indepen- 
dent quantitative classification. In this way, we can determine which rilles 
are considered by most workers to be sinuous and determine their structural 
features quantitatively. In addition, the quantitative results might be used 
to further differentiate sinuous rilles into other subclasses. 
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FINITE FOURIER APPROXIMATIONS OF LUNAR RILLES 


If one end of a lunar rille is connected to the other by a straight line, 
this base line can be considered the X axis of a Cartesian coordinate system. 
The course of either bank of the rille can be considered mathematically as 
F(X). In thi s case, F(x) is usually an irregular single-valued function over 
most of the rille length. Most single-valued functions can be represented by 
the Fourier series if an adequate number of harmonics is considered. However, 
rilles often recurve upon themselves and in this interval of X there is more 
than one value of F(x) and the simple Fourier series does not define the. func- 
tion. In general, this interval is small compared to total rille length so 
that the discontinuity has small effect on the spectral components of the rille 
and the method is suitable for defining a rille quantitatively. 


Data for analysis of a rille is obtained by measuring discrete values of 
F(X r ) at N+l equally spaced points along the base line, l, of a given rille. 
Thus, the rille is represented by the data function F(x r ) where X r = r£/N, 
r = 0,1,. . . , N and F(0) = 0. This data function can be approximated by any 
given set of harmonics n < N at each of the N+l points by the finite 
Fourier series : 

n sin X kir 

F n< x r) = l \ T— i » = 1.2,. . . ,N 

k=l 

The Fourier coefficients, B^, are determined by least-squares estimates: 

0 N sin X kiT ' 

A — I F(X r> f~ 

r=l 


When n = N the Fourier approximation, Fjj(X r ), is equal to the empirical 


representation of the rille at each of the data points, 
in fit is estimated by: 

N [F„(X ) - F_(X T .)] 2 
E(n) = l 
x=o 


'N' 


N + l 


When 1 < n < N error 


The method used to estimate the least number of harmonics n that provide an 
adequate fit was to compare the percent reduction in error resulting from an 
increase in the number of harmonics with the reduction in error that resulted 
from a preceding increase in the number of harmonics; that is, comparisons of 
the type: 

[E(n)-E(n+A)/E(n)] * 100 


were computed for n = A, 2A , 3A, . . . until the percent reduction was less 
than 1. These calculations were also analyzed with regard to the average 
amplitude of each rille. 


EXAMPLES OF SPECTRAL COMPONENTS 0^ LUNAR RILLES 


The method of Fourier analysis of rilles is ideally suited for an objec- 
tive classification of lunar rilles because it can define quantitatively all 
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planimetric features of rille shape except recurvature. A modification of 
the simple Fourier analysis is being developed to allow for the description of 
recurvature. The degree of approximation of a given rille by different sets 
' of harmonics of the Fourier series is shown graphically in figure 1. Approxi- 

mations to the rille at the top of the figure result from the fewest number of 
harmonics. As higher order harmonics are included, the calculated values of 
F(X) more nearly approximate rille 22 of the catalog shown at the bottom of 
the figure. Inspection of figure 1 shows that approximation of the rille 
using the first 109 harmonics is very accurate. 

There is a wide range in pattern of harmonics associated with lunar 
rilles . It is this variation of values of Fourier coefficients and variation 
in significance of harmonic order that suggests that rilles can be classified 
using this technique. Figure 2 shows the amplitudes of % of the three 
rilles (rilles 4, 12 , and 22 of the rille catalog). Amplitudes of are 

plotted as a function of harmonic order. The wavelength of the first harmonic 
is twice the rille baseline length, the wavelength of the second harmonic is 
equal to the baseline length, and so on. The number of harmonics with signifi- 
cant amplitudes that are required to define each of the rilles is quite dif- 
ferent, as are the amplitudes. While rilles 4 and 22 are characterized by a 
small number of harmonics, rille 12 is adequately described only by a complex 
spectral pattern that indicates the importance of higher order harmonics. The 
wide range in patterns of spectral components indicates that lunar rilles can 
be easily classified using Fourier analysis. 


RILLE CATALOG 


Lunar Orbiter photographs were examined for rilles shown at resolutions 
suitable for analysis. Most of the photographs in the catalog are taken from 
LO IV, although some are from L0 V. The rille catalog is in two parts. The 
first section (rilles l-9l) shows simple rilles (those that usually contain no 
tributaries; however, some contain one or two branches or tributaries). They 
are obviously isolated rilles that are not part of a rille system. The second 
section (rilles 1001-1029) shows complex rille systems. These sometimes occur 
completely within craters and they are characterized by many intersecting or 
branching sections of rilles. 

Each plate of the catalog shows one or more rilles or rille systems and 
the latitude and longitude of the center of each. Each rille is numbered, 
rille branches are lettered, named rilles are labelled, and a scale is given 
on each photograph. 

Although the catalog was produced primarily as a working document for a 
quantitative analysis of lunar rilles, it may also be useful as a reference 
for investigators making other quantitative and qualitative studies of rilles. 
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LO V - M-190 
43° W, 27° N 
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HADLEY RILLE 
LO IV - H-102 
3° E, 25° N 
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RIMA CONON (11) 
LO IV - H-102 
2° E, 19° N 
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RIMA BOND I 
LO IV - H-79 
36° E, 32° N 
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ALPINE VALLEY 
LO IV - H- 115 
3° E, 49° N 
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LO IV - H-183 
79° W, 34° N 
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15 km 



LO IV - H-78 
29° W, 15° N 
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RIMA THEAETETUS 
LO IV - H-103 
6° E, 33° N 



LO V - M-130 
4° W, 49° N 
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LO IV - H-133 
33° W, 8° N 
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LO V - M-166 
31° W, 13° N 
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LO IV - H-151 
47° W, 29° N 
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LAVOISIER CRATER 
LO IV - H-189 
81° W, 38° N 
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PITATUS CRATER 
LO IV - H-119 
14° W, 30° S 
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MERSENIUS CRATER 
LO IV - H-149 
49° W, 21° S 
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ATLAS CRATER 
LO IV - H-79 
45° E, 47°N 








RIMA SIRSALIS (1016) 

LO IV - H-156, 161, 168 
65° W, 12° S 
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100 km 
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ALPHONSUS CRATER 
LO IV - H-108 
3° W, 13° S 
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HUMBOLDT CRATER 
LO IV - H-27 
80° E, 28° S 
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n /X r kTi\ 

F n (X r ) = l B k sin l-±— ); n = 1,2,. . .,N 

k=l V ' 

Page 2, second equation should read: 


B k 


2 

n 


N 

l 

r=i 
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